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We study linear perturbations around static, spherically-symmetric spacetimes in f{R, C) gravi- 
tational theories whose Lagrangians depend on the Ricci scalar R and the parity violating Chern- 
Simons term C. By an explicit construction, we show that the Hamiltonian for the perturbation 
variables is not bounded from below, suggesting that such a background spacetime is unstable 
against perturbations. This gives a strong limit on a phenomenological gravitational model which 
violates parity. We also show that either R = const or g^^ = is a necessary and sufficient condi- 
tion for the stability. We then implement in detail the perturbation analysis for such theories which 
satisfy the stability conditions. For (.>2, where l is the usual integer for the multipole expansion, 
the number of propagating modes is three, one from the odd and the other two from the even, all of 

O ' which propagate at the speed of light. Unlike in the case of f{R) theories, these modes are coupled 

yJ , to each other, which can be used as a distinctive feature to test the parity violating theories from 

II ■ observations. The no-ghost conditions and no-tachyon conditions are the same as those in f(R) 

theories. For the dipole perturbations, the odd and the even modes completely decouple. The odd 
mode gives a slowly rotating black hole solution whose metric is linearized in its angular momentum. 
We provide an integral expression of this solution. On the other hand, the even mode propagates 
at the speed of light. For the monopole perturbation, in addition to a mode which simply shifts the 

O \ mass of the background black hole, there also exists one even mode that propagates at the speed of 

O^i light. 

iH . 

I. INTRODUCTION 

^ . . . ^ 

^ ' General relativity (GR) has been frequently tested both experimentally and observationally over many decades to 

lO \ check whether it is really the correct theory of gravity [l|. In the weak gravitational field regime, like in the solar 

^^ ■ system or on the Earth, no deviations from GR have been detected so far. In this sense, GR is a successful theory 

^^ of gravity. However, the need for the introduction of dark energy to explain accelerated expansion of the universe 

■ (for recent reviews on dark energy, see e.g., Refs. @, 9) ™^^y be a signal that GR breaks down in a regime where 

C ■ ' the strong field effects become important. Furthermore, although there are few observational tests of GR in the 

^^ , strong field regime, in the future, we will be able to probe the nature of strong gravity, for example, by observing 

, - the gravitational waves coming from the vicinity of black holes (BHs) [J, [a| . These facts have provoked alternative 

. . theories of gravity such as f{R) theories (for recent reviews of f{R) theories, see e.g. Refs. 0, [tI) and have led us to 

^ understand theoretically what kinds of different phenomena are expected in such theories [8l-[l9| . 

In light of this situation, it is interesting to consider gravitational theories which violate parity due to the so-called 
Chern-Simons(CS) term, or the Pontryagin density, C = ^eap-ysR"'^ i^^^R'^^^'^ , where eap-ys is the totally antisymmetric 
cd ' tensor. The presence of the t^p-ys tensor manifests the parity violation. Because the CS term can be expressed as a 
divergence, simply adding C into the Einstein-Hilbert action amounts to an addition of a total derivative term in the 
action and does not change the theory. A gravitational theory in which the CS term is coupled to an external scalar 
function was introduced in Ref. [20| • The external scalar function was promoted to a dynamical field in Ref. [21| and 
many studies based on such theories have appeared [22h35I |. For a recent review on the Chern-Simons gravity, see 
Ref. [36] . Noticeable features of the CS term are that it identically vanishes for the Friedmann-Lemaitre- Robertson- 
Walker (FLRW) metric and even for the scalar type perturbations on top of it, and for the spherically symmetric 
metric. As a result, the cosmological and the Solar System constraints achieved so far constrain the CS gravity only 
very loosely. Therefore, we need to go beyond such simple spacetimes to make the lurking CS term yield different 
phenomena from GR. For example, the CS term enters the game in the tensor perturbations on the FLRW background, 
which suggests that observing gravitational waves is an effective approach to test the CS-gravity ,37.-:4S] . Another 
useful approach is to consider perturbations around the static, spherically-symmetric background such as Schwarzshild 
BH, where the CS effects are expected to show up [22 - [25l . [2a - l30l . |33] . 

In this paper, we consider the gravitational theories whose Lagrangian is a general function of R and C, f{R, C), 
and develop linear perturbation theory around the static and spherically symmetric spacetime. Similar analysis has 
been done for f{R,Q) {Q is the Gauss-Bonnet term) theories in Ref. [44]. Unlike in the case of f{R) theories which 
can be mapped into equivalent theories where a scalar field having self-interacting potential is minimally coupled to 



X 



Einstein-Hilbert gravity, f{R, C) theories cannot be mapped into theories where C is coupled solely to a dynamical 
scalar field due to nontrivial transformation property of the CS term under the conformal transformation. Our aim is 
to clarify both quantitative and qualitative behaviors of the perturbations. To be more precise, we will derive no-ghost 
and no-tachyon conditions which are necessary to ensure stability of the background spacetime against perturbation, 
obtain dispersion relations for the propagating modes and find features that are characteristic to parity violating 
theories. 

The BH perturbation for the so-called nondynamical Chern-Simons theories where the Lagrangian consists of the 
Einstein-Hilbert term plus the CS term C multiplied by a nondynamical scalar field has been studied in Ref. [33|- 
As we will show in the next section, the general f{R, C) theories we consider do not cover such nondynamical CS 
theories. Correspondingly, the results presented in Ref. [33[ cannot be applied to our study. 

The organization of this paper is as follows. In Sec. |lTl we derive the gravitational field equations and apply them 
to the static and spherically symmetric background. In Sec. IIIIl we will develop linear perturbation analysis on that 
background by expanding the action in second order perturbations and show that the perturbations are unstable in 
general. Then in Sec. IIV| we will study in detail a special class of theories in which the instability of the perturbations 
is avoided. The last section is the conclusion. 

II. BACKGROUND EQUATIONS 

We study f{R, C) theory, where the action is described by a general function of Ricci scalar R and CS term 

C = ^eap-fSR°'^fj,i,R'^ ^" , 
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Here, Mp = 1/^/EttGn ~ 4.34 x 10 ^g is the reduced Planck mass. We can rewrite the action ([T]) as 



A/2 , 



S ^ ^ J d^x V^ iRF{X, s) + WiX, s)C - V{X, s)) , (2) 

where A and s are auxiliary fields and 

F{X,s) = ^^^, W{X,s) = ^^^^, V{X,s)^XF{X,s)+sW{X,s)-f{X,s). (3) 

The variations with respect to A and s yield the following constraints: 

X^R, s^C. (4) 

By substituting them, it can be confirmed that the second action ([2]) actually reduces to the first one ([!]). The 
gravitational field equations obtained from the action ^ are given by 



1 1 



V^V.F - g^^DF - 2e''\ (R%)afi^.^-fW + 2V^W V„i?,)^) - ig^.y 



0. (5) 



Here, the parentheses around the indices indicate symmetrization. Si ,,,a = {S^^ + S^^)/2. 

Let us show that the nondynamical CS theories considered in Ref. [33|, whose action is given by 



^non-dynamical ^^ ^^^ V^ {R + WC) , (6) 

are not covered by Eq. ^. In Eq. ([5]), VF is treated as a field. Apparently, we can reproduce Eq. © by setting F = 1 
and y = in Eq. ([2]). However, by using the definitions of F, W and V given by Eqs. ([3]), the conditions F = 1 
and V = lead to VF as a constant. Hence, the action describes GR instead of Eq. ([6|). Thus, the nondynamical CS 
theories are not included in f(R,C) theory, and vice versa. Clearly, dynamical CS theory, whose Lagrangian has the 
kinetic and potential terms in addition to Eq. ([5]), is also not included. The analysis for nondynamical and dynamical 
CS theory is given in Ref. [45| . 

Throughout the paper, we consider a static and spherically symmetric spacetime as a background, i.e., we set the 
background metric as 

ds^ = -A(r)dt^ + —— + r^ (de^ + sin^ 9 dip^) . (7) 



For this metric, nonvanishing and independent elements of X^^ are given by 

A{r{rB'F' + 2rBF" + ABF' + rV) + 2F{rB' + B -I)) 

^" - ^^r^ ' («) 

rBA' [rF' + 2F) + A{r{ABF' + rV) + 2{B-l)F) 
''^ ~ 2t-'ABF ' ^^^ 

r {-rBFA!"^ + A {rFA'B' + 2B {rA'F' + F {rA" + A'))) + 2A^ {rB'F' + FB' + 2rBF" + 2BF' + rV)) 
Xee = ^^f ^10) 

where ' means derivative with respect to r. The function W does not appear in X^^ because of the spherical symmetry 
of the background. We can convert the background equations Xtt = 0, Xrr — and Xgg = to express A" , F" and 
V in terms of A, B and F and their first derivatives, 

,_ -rAA' {rFB' + 2B{rF' +F))+r^BFA'^ + 2A^{F{rB' + 2B-2) + 2rBF') 
^ ~ 2r^ABF ' ^^^^ 



„ _ {rF' + 2F) jBA' - AB') 
2rAB 



(12) 



rBA' {rF' + 2F) + 2A {2rBF' + {B - 1)F) 

We use these equations when the second derivative terms appear as a result of integration by parts to rearrange terms 
in the second order action for the perturbation. 

III. PERTURBATION 

In this section, we calculate the perturbative action. We start from the metric perturbation with a brief review of 
the Regge- Wheeler formalism. As a result of the parity violation term, the odd and even modes do not decouple from 
each other within the system of equations and we have to deal with it all together. We derive the second order action 
and find that there exists an instability due to the presence of the higher derivative terms. We assume i > 1 when 
we construct our final action after the elimination of auxiliary fields. There is also an issue regarding residual gauge 
degrees of freedom for ^ = 0, 1 and we will consider those modes separately in the next section. 

A. Decomposition of the metric perturbation 

Before studying the metric perturbation of a spherically symmetric static s pac etime for f{R, C) theories, let us 
briefly review the formalism developed by Regge and Wheeler [4^, and Zerilli [43] to decompose the metric pertur- 
bations according to their transformation properties under two-dimensional rotations. Although Regge, Wheeler and 
Zerilli considered the perturbation of the Schwarzschild spacetime (namely in GR), the formalism solely relies on the 
properties of spherical symmetry and can be applied to f{R, C) theories as well. 

Let us denote the metric, slightly perturbed from a spherically symmetric static spacetime, by ^^j/ = (7° -I- /i^j/. 
Hence hf^^, represents infinitesimal quantities. Then, under two-dimensional rotations on a sphere, htt,htr and hrr 
transform as scalars, hta and hra transform as vectors and hab transforms as a tensor (subscripts a, h are either 9 or 
if). First, any scalar s can be decomposed into the sum of spherical harmonics as 

s(t, r, e,^)=Y^ Sim{t, r)Yi^{0, cp). (14) 

£,171 

Second, any vector Va can be decomposed into a divergence part and a divergence-free part as follows: 

K(t,r,0,^)=V,$i+S^V6$2, (15) 



where $1 and $2 are scalars and Eab = vdet7 eat with 'jab being the two-dimensional metric on the sphere and Cob 
being the totally antisymmetric symbol with ee^ = 1. Here Vq represents the covariant derivative with respect to the 
metric "fab- Since Va is a two-component vector, it is completely specified by the potentials $1 and $2. We can then 
apply the scalar decomposition (J14p to $1 and $2 to decompose the vector quantity Va into spherical harmonics. 
Finally, any symmetric tensor Tab can be decomposed as 

Tab{t, r, 0, V3) = VaVb^i + 7a(,*2 + ^ (^^a^VeVfe^g + £;fc"VeVa*3) , (16) 



where \l/i, ^2 and ^3 are scalars. Since Tab has three independent components, ^1, ^^2 and ^3 completely specify 
Tab- Then we can again apply the scalar decomposition p^ to ^1, ^2 and ^3 to decompose the tensor quantity into 
spherical harmonics. We refer to the variables accompanied by Eab as odd-type variables and the others as even-type 
variables. 

Let us now apply this decomposition to the metric perturbation /i^^ = 5^1, — gliJ . For the odd-type perturbations, 
they can be written as 

htt^O, htr = 0, hrr^O, (17) 

hta = J2 '^oMt, r)Eabd^Yim{e, ¥>), (18) 

hra = Y. h^-imit, r)Eabd''Yi„,{e , (^), (19) 

e,m 

h-b^lj2 h2,lmit. r) [Ea'V,WbYt^{e, ip) + E^^'V ,W aY,m{e , 'p)] ■ (20) 

For the even-type perturbations, they can be written as 

htt = A{r)Y,HoAm{t,r)Yi,,,{e,if), (21) 

htr = Y. Hl,tra{t, r)Yt„,{e , If), (22) 



e,m 



--^y^H2.Ut.r)Y,m{e,v). (23) 

hta =Y,Pi-^{t,r)daYira{0,^), (24) 

Ka = ^aim{t,r)daYtm{0,(p), (25) 

hab = J2 ^^'"(*' r)9abYe,nid, ^) + J^ '^''™(^' ^)'^ a^ bYemid , v) ■ (26) 

Because of general covariance, not all the metric perturbations are physical in the sense that some of them can be set 
to vanish by using the gauge transformation x^ — >■ x^ -1-^'', where £,^ is infinitesimal function. A gauge transformation 
acting on the odd-type perturbations can be written as 

^a = Yl ^^"(^' r)Ej'dbYirn{0, v) , (27) 

where Afm(i,r) are arbitrary functions. Gauge transformations on even- type perturbations can be written as 

6 = 5I^^'"(*'^)^^™(^''^)' ^'■-^i?£r„(i,r)y,™(^,(^), ea = ^e,™(t,r)aar,™(^,(^), (28) 

£,m £,m £^m 

where Tim{t,r), Tim{t,r) and Qim{t,r) are arbitrary functions. Under these gauge transformations, we can find the 
transformation rule for each metric component. For the odd-type perturbations, we find 



ho,ern{t, r) -> hojrnit, r) + Aim{t, r), (29) 

2 
hi,im{t, r) -^ hi,e,n{t, r) + A^„(t, r) - -Aim{t, r), (30) 

r 

h2,emit, r) ->■ h2,imit, r) + 2A<>,„(t, r), (31) 



where a dot denotes the time derivative. For the even-type perturbations, we find 

2 • A' B 

Ho,emit,r) -^ Ho,imit,r) + —Tim{t,r) — i?^„j(i,r), (32) 

A' 

Hijyn{t,r) -^ Hi^em{t,r) + Remit, r) +Tl>„^{t,r) - —Ti„i{t,r), (33) 

H2Mt,r) ^ H2,im{t,r) + 2BR'e,Jt,r) + B'Re„r{t,r), (34) 

I3em{t, r) -^ Ptm{t, r) + T^„(t, r) + Qi^(t, r), (35) 

2 

(^£m{t,r) -^ ae„i{t,r) + Ri„i{t,r) + Q'f„^{t,r) efo„(t, r), (36) 

r 
9 R 

Kemit,r) ^ Kt,,nit,r) + — Remit, r), (37) 

r 

Gim{t, r) -^ Gernit, r) + 2e^„(t, r). (38) 

Since no derivative appears in the transformation rule for /i2 , we can completely fix a gauge by imposing a condition 
/i2 — 0. This gauge fixing is called Regge- Wheeler gauge. For the even-type perturbations, complete gauge fixing is 
achieved by imposing ,9 = 0, K = and G = 0. We will use these gauge conditions in the calculation of the second 
order action. 

We remark here that the parity violation induces the coupling between odd and the even modes. This can be 
understood by comparing f(R, C) gravity with its special case, f(R) gravity. In f{R) gravity, where / depends only 
on R and the theory does not violate parity, the linearized equations of motion (or equivalently, the second order 
action) for h^^ can be decomposed into ones that only contain odd-type perturbations and ones that only contain 
even-type ones. This decoupling drastically simplifies the perturbation analysis since it allows us to consider equations 
(or an action) with fewer perturbation variables. However, this decomposition does not hold for f{R, C) gravity due 
to the explicit violation of parity. Therefore, we must treat both the odd and the even modes at the same time, which 
we will do in the following sections. 

B. Second order action 

In addition to the metric perturbations, we also need to perturb the other functions A and s that appear in the 
action. Just for later convenience, instead of perturbing A and s as the fundamental fields, we treat 5F and 5W as 
perturbation variables. These fields must be also decomposed into the spherical harmonics, 

6F^Y.^^^ra{t,r)Y,m{e,ip), (5VF = ^ JVF,™ (t, ^Ff™ (0, (^). (39) 

The relation between {5F, 5W) and (5A, 5s) is given by 

5F ^Fx5X + F,5s, 5W ^FM + W^Ss, (40) 

where we have used the identity Fs = Wx = Qxds *° replace W\ with Fs in the second equation. 

With these perturbation variables, expanding the action ^ to second order (the first order part automatically 
vanishes because of the background equations) yields the following action. 



/■ 



dt dr C, (41) 

where C is written as 

C = Ho{aiH2 + 02^2 + "3^ + a4a' + a^SF + aeSF' + ajSF" + a^ho + agh'^ + aio/iQ + aii^i + ai2h[) 
+biHl + Hi{b25F' + h^5F + biH2 + ha + hliQ + b-jh'Q + bg^i + ^g/^i) 
+ciHl + H2{c2a + C'iSF' + c^6F + c^h'^ + c^h^ + cjhi) + csH2S'F 
+dia'^ + d20? -t- a{dy,5F' + d^SF + d^h'^ + dg/io + d-jh'^ + dgho + dghi,ttt + rfio^i 
+eihl + ho{e2hi + esSW' + e^SW) + e^h'^hi + e^h'^ + ejhl + eg/ii + eghiSW 
+fidF^ + f2SF6W + fsSW^ . (42) 



Since different {I, m) modes do not mix with each other, we pick up particular {£, m) modes. Because of the spherical 
symmetry of the background spacetime, the action for m ^ Q modes takes exactly the same form as that for m — Q, 
which enables us to set m = without loss of generality. From now on, we abbreviate the subscripts £ and to. Explicit 
expressions for the background-dependent coefhcients are given in the appendix [X] We first notice that, due to the 
parity violating nature of the CS term, there appear mixing terms of odd and even perturbations in the action. For 
the mixing terms that do not contain 5W , all the coefficients vanish if W{r) does not depend on r. This reflects the 
fact that the Chern-Simons term C is a total derivative and affects the field equations of motion only when W , which 
multiplies C, depends on r. 

This action shows that not all of the variables are dynamical. Actually, we see that Hq, Hi and 5W are auxiliary 
fields. Therefore, they can be eliminated from the action by using their equations of motion. Since iJg appears only 
linearly, the variation with respect to it gives a constraint among the other fields. We eliminate H2 by using this 
constraint. The variation with respect to Hi and SW gives the equation of motion of each variable, respectively. After 
substituting the constraints and many integration by parts, we end up with the following Lagrangian density: 

C ^ pi'hl +P2hi{rh'Q - 2ho)+p3h'^ + Pihl + pslij + peSF^ +pj$'^ + psho6'F + pglioP + pio$6'F + pnh'^ 
+Pi26F'^ +Pi3^'^ +Pi4h'oSF' +pi5h'ol3' +piel3'SF' +pijh'ohi +pi8hohi + p ig h'^S F +p2oh'o (3 

+P2lhidF + P22hlP + P23SFI3' + P2ihl + hn{p25SF + P2&P) + P2Thl + P28SF^ + P29SFP + P30/32. (43) 

Since all the fields (/iq, /ii, /?, 6F) have time derivatives that are not removed by any integration by parts, all of 
them are dynamical fields. Hence, this is our final Lagrangian. Explicit expressions of the background dependent 
coefficients pi,- ■ • are given in the appendix |51 However, since most of pi, • • • have very long expressions, we only 
provide them in a form truncated at the leading order in F', W, Fx, Fg and Ws all of which vanish in the GR limit. 
In addition to the mixing terms of odd and even modes, there is another big qualitative difference between the 
above Lagrangian and that in f{R) gravity. The above Lagrangian contains a term h^. This term results in fourth 
order differential equations for hi with respect to time, 

'2pi-g^+P2 (^-5^-2-^1 ~ 2p5hi - pnh'o + pisiio - P2iS'F - p22^ + 2p27hi =0. (44) 

In f{R) gravity, there are no second derivative terms in the Lagrangian and the resulting equations of motion are 
second order. We can confirm this fact by looking at the explicit expressions for pi, p2 and pa, which are given by 

327r£(£+l)M|,iy'2 2pi 
Pi = 777175"' P2 = , P3=Pi- 45 

For f{R) gravity, we have W = 0. Therefore, pi, p2 and p^ vanish and the Lagrangian does not contain second 
derivative terms. 

The presence of the hi term in the general f{R,C) gravity is a signal that the theory is plagued by instability 
\4§\ . Actually, we can show that Hamiltonian corresponding to Eq. (|43|l is not bounded from below. To construct the 
Hamiltonian, we find it useful to introduce a new field q and rewrite Eq. (j43]) as 



£ = -Pi 



q'^ + 2q(hi-h'g + -ho 



oiPiV 4pi 



hl+p5hl+pe6F + prP^ + pghoSF + pghoP 



+Pwp5F + piih'i^ + pi2(5F'2 + pi3/3'^ + Puh'oSF' + pi^h'^p' + pieP'SF' + pnh'Jn + pishohi + pigh'gSF 
+P2ohol3 + hi{p2i5F + P22P) + P23SFf3' + p2ihl + ho{p25SF + p2e/3) + P27h\ + P2?.5F'^ + p29<5i^/3 + P3o/?t46) 

where we have used relations ps)) . It can be confirmed that Eq. (pS)) reduces to Eq. (H5)) after eliminating q by using 
its equation of motion. The new Lagrangian (j46l) is much more familiar than Eq. (1431) since it does not contain either 
the second time derivative term nor the mixing derivative term H'q. Now, we can use the standard canonical formalism 



to construct the Hamiltonian. The conjugate momenta are defined by 



^ = -2p,(^h,-h', + ho], (47) 



P4 - 2 — ' 
r 



ho+p8SF+pgf3+pishi, (48) 



oho 

dC 
TTi = -^ = -2pig + 2p5hi + pirh'o + p2idF + p22/3, (49) 

ohi 

dC ... 

TTF = --^ = 'ipeSF + p^ho + pio/3, (50) 

9£ • ■ 

TTi3 :^ —J- =2p7/3+pgho+pioSF. (51) 

op 

We can always solve these equations in terms of the field time derivatives since a Jacobian of these transformations 
is given by 



2 



det ^. „! I = 5-, (52) 

\dq^dqjj {2i+lfA^F {rBA' {rF' + 2F) + 2A{F {-~2B + l"^ + i) - rBF')y 

which is not zero in general (the special case where VK is a constant will be considered later). Here we defined 
<7i = q, q2 = hf), 53 = hi, qi = 5F, (75 = /3. Therefore, the Lagrangian is not singular and there are no primary 
constraints among the canonical variables. The Hamiltonian is then given by 

H = dr v. = dr [TTqq + ttq/io + tti/ii + iTpSF + -Kfjfi — C\ . (53) 

It can be confirmed that a matrix Kij , defined by the momentum part of the Hamiltonian, 

n D K.jTT.TTj +■■■ , (54) 

where • • • represents terms that are not quadratic in tt^, has a vanishing component for K33 and a nonvanishing 
component for Ki^, 

K33 = 0, i^i3 = --L. (55) 

4pi 

In addition, we find {tti, tt^} and {112, 7^4, tts} are decoupled from each other. Combined with these facts, the 
subspace of the Hamiltonian spanned by tti and tt^ yields a negative determinant of the corresponding subkinetic 
matrix, 

KnKss - Kf, = -Kf, < 0. (56) 

This means that the Hamiltonian can take arbitrary negative values by suitably choosing the values of tti and 773 and 
hence it is not bounded from below. This result shows that the general f{R, C) gravity has the problem of having a 
ghost around the static and spherically symmetric background and provides a severe condition on the functional form 
of/(i?,C). 

There are two possible cases where the presence of the ghost does not become problematic. The first one is to 
assume that the f{R, C) theory under consideration is an effective theory which is valid only on length scales larger 
than a certain length dc- From this point of view, the presence of the ghost does not matter if its mass is larger than 
the energy scale d~^ since the dynamics of the ghost cannot be described by the low energy /(_R, C) theory. The 
more fundamental theory which is valid above d~^ may cure the problem. The mass of the ghost can be evaluated 
as follows. Neglecting the gradient terms for the fields which are not important here, the Hamiltonian density can be 
written as 

% = KijTTiTTj + BijTTiqj + Mijqiqj. (57) 

There are mixing terms between tt^ and qj that are represented by a matrix Bij. Let us first eliminate the mixing 
terms by the following canonical transformation: 

^T,=P^-l{K-'B)^^Q„ (58) 

q^ = Q^. (59) 



In terms of the new canonical variables, the Hamiltonian density can be written as 



n = K.jPiPj + M - -B'^K-^B ) Q,Qj. (60) 



1 

4' 

It turns out that the new mass matrix does not mix {Qi, Q3} and {Q2, Q4, Qs}- Therefore, the first and the third 
canonical fields form a closed system. The corresponding sub-Hamiltonian is given by 

Hsub = Ki^Pi + 2Ki3PiP3 + MiiQl + M33QI (61) 

where M = Af — jB^ K^^B and each clement is given by 

{2e+l){6M{e+l)F'^F^-r*Fx) 2£ + 1 



20487r^(£+l)Af|,r4i^AW^"' ' 1287r^ (^ + 1) M|.Vl/'2 ' 

" 2^+1 ' '^ 32 (2^3 + 3^2 _ 3^ _ 2) r8F2l^'2 ' ^ ' 

We notice that both Mn and M33 are negative definite, which means both Qi and Q3 become tachyonic as well 
if the magnitudes of their masses are smaller than d~^. We can make the sub-Hamiltonian, which is a sum of two 
independent harmonic oscillators, by the following canonical transformation. 



Pi = \/-Mii (Pi cos6 + P3 sinS) , P3 = V-M33 {-Pi sinS + P3 cos (5) , 

Qi = ^^^^ ((5iCOS(5 + (53sin(5) , (53 = ^^^^ (-Qi sin(5 + (53 cos(5) , (63) 



where 6 is determined from the equation 



2X13 /M33_ / 3 _ ,4^ ^^,,^ , ^^^2W2^ 



tan 2^ --^^J^^- V £^T7^ ^^^('^'^^ - 64£(^+ DFf^-) . (64) 

Here, the sign function is defined as sgn(a:;) = +1, 0,-1 for x > 0, a; = 0, x < 0, respectively. The new sub-Hamiltonian 
is then given by 



Hsuh ^ ~ cosS (KnMiicosS + 2Kl3^/MnM33smS) P^ -sinS [KnMiismS - 2Ki3y/MnM33COsS) P^ ~Ql-Qi 



The coefficients in front of Pi and P3 must be larger than d^ ^ jn order for those fields to be in the high energy regime 
where the effective f{R, C) theory does not work. This leads to a condition. 



\KiiM: 
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^Fx-6Uii+l)F!F- ^^_,^ ^^^^ 



QAr^FxW^ 

In particular, when the second term in the numerator is negligible, we obtain the very simple condition for W' , 

\W'\ < 4. (66) 

Since W', which has dimensions of length, represents how large the effects of the Chern-Simons term arc, this condition 
says these effects are suppressed on distances larger than dc- 

The second possibility where the presence of the ghost does not become problematic is that f(R, C) belongs to the 
special class in which W' = is satisfied identically. Using the background metric, W' can be written as 

W = FsR' + WsC = FsR', (67) 

where we have used an identity C — for the background metric. Therefore, if f{R,C) satisfies either P^ = or 
R — const, we have W' — identically. For example, P, = is trivially satisfied if f{R, C) takes a separable form, 
i.e., f{R,C) — fi{R) + /2(C), where /i and /2 are arbitrary functions of R and C, respectively. f{R) gravity is 
included in this case. The second case R = const is satisfied, for example, if the Schwarzschild metric is a solution of 
the model. In either case, we have still many f{R, C) theories. We deal with this class of f{R, C) theories in the next 
section. 



IV. STUDY OF SPECIAL CASES WITH W = 0. 

As we have shown in the previous section, the general f{R, C) theories with nonvanishing W' have the problem of 
instability. Thus, the cases with W = are more phenomenologically interesting and deserve further investigation. 
In this section, we study the second order perturbation again for this case. We clarify the number of propagating 
modes and derive dispersion relations for them. First, we analyze the general case for ^ > 2, and consider the special 
case for £ = and 1=1. 

A. Second order action again 

The Lagrangian for the general /(-R, C) theories (j43|) can be also used for the special case W = 0. Since the 
condition W' = makes some terms identically vanish and the resulting Lagrangian is greatly simplified, it is better 
to write the simplified Lagrangian, 

C = qi [K -hiY + q25'F^ + qzSF'^ + QiSF/S + q^SF'P' + 96/3' + qjP'^ + q^hoh^ 

+q9SF{h'o - hi) + qw(35F' + quh^ + quh^SF + qishj + qi^SF^ + qi^(35F + qieP\ (68) 

where the background-dependent coefficient can be read from pi, ■ ■ ■ by imposing a condition W' = 0. The odd and 
even modes are still coupled. The background-dependent coefficients responsible for the coupling are given by 

IG-Kiii + l)MlFs (rBA' {rF' + 2F) - 2A {rBF' + 2{B - 1)F)) 2 

*- i2i+l)r^AFF, ' '^12 ---99. (69) 

Interestingly, the odd and even modes are decoupled if Fg = 0. Furthermore, if Fg = 0, we verified that all the 
coefficients qi for the even modes are exactly the same as those for the f{R) theories. On the other hand, the 
coefficients g^ for the odd modes depend on Wg- Therefore, for the f{R, C) theories having the property Fg = 0, the 
signature of the parity violation appears only in the odd modes. In the following analysis, we do not assume Fs = 
to keep the procedure as general as possible. 

Similar to what we did in the previous section, we can introduce a new variable q to write the above Lagrangian as 



£ = qi 



2q[h'o-hi + -ho ] - q^ 



q2S'F^ + qsSF'^ + q^SF^ + q^SF' 13' + qefi^ + qr/S'^ 



+qsSF(h'o - /li) + qiofiSF' + (qn - ^q's - ^) h^ + qi2hoSF + qishj + q^SF^ + q^^pSF + q^eP^. (70) 

Here we have used a relation qs = 2qi. We can verify that the new Lagrangian reduces to the original one by 
eliminating q using its equation of motion. The point of the new Lagrangian is that it contains derivatives of ho and 
hi up to at most first order while the original one contains h'^^ and hf. Therefore, by doing integration by parts, we 
can rewrite the new Lagrangian in such a way that any derivative of ho and hi does not appear anymore. After this 
procedure, both ho and hi become auxiliary variables. For i > 2, quadratic terms in ho and hi exist and we can 
eliminate ho and hi by using their equations of motion. After this, we finally obtain a Lagrangian which contains 
only q, SF and f3. We find that the final Lagrangian can be formally written as 

C = kijqiqj - dijq^q'j - Ci^q'^qj - mijqiq.j, (71) 

where we have defined (qi, ^2, 93) = {5F, /3, q). The cases where ^ is either or 1 will be studied later. 
A determinant of the kinetic matrix kij is found to be 



det(fc.,) = .,,3/2 "^ J' (72) 



(2£ + IfA^F^ {^y^ {rBA' (rF' + 2F) + 2A (F (-2B + P + £) - rBF')) 
where Y is defined by 



mije + 1) (F^Ws - F^) {rBA' {rF' + 2F) - 2A {rBF' + 2{B - l)F))^ ^ ^, ^3 
^= ^T^r +AF. (73) 

This is not zero in general. Therefore, all the variables are dynamical and there are three propagating modes, one of 
which is odd (i.e., q) and the remaining two are even {i.e., SF and /?). This structure is the same as that of the f{R) 
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gravity theories where there is one propagating odd mode and two propagating even modes. This result shows that 
the condition W' — kiUs all the pathological modes which, as we found in the previous section, exist in the general 
f{R, C) theories. To see if the modes are ghosts or not, we find it convenient to evaluate ^33 and k22k33 — fc^g, 



«33 



27ri{i + l)2Af|,r2y2 



(2£4 + 5^3_5£_2)^5F^ 

32t:^£^{£ + if MJ^r'^Y'^ B"^ {{2^ + 2i~l) r^F'^ + 2{e+l~2) rFF' + 2{e+i-2)F'^ 
{21 + 1)2 (£3 + 2^2 _ £ _ 2) ^4i^4 (j,_g^/ (^p, _^ 2F) + 2A {F {-2B + P + () - rBF')f 



(74) 



To avoid the ghost, we need to impose F > 0. ^ With this condition, both dei{kij) and fcgg are clearly positive 
definite. After a short calculation, it can be verified that ^22^33 — ^23 is also positive definite. Therefore, as is the 
case with f{R) theories, F > is the no-ghost condition for /(i?, C) theories that satisfy W — 0. 
We can derive the dispersion relations for the three modes from an equation, 

dei{-u'^hj + k'^dy) = 0. (76) 

Explicit calculation gives 

7Q8nH^ i-e -2e^ + e + 2? M^rBF^qf.f^ (uo^ - PABf 
det(-o.2^ ^ ^2^ ) ^ ^ '_ ^_±^ '_ (77) 

{2(. + lYqi3{2rqii-rq'^-2qs){rBA' {rF' + 2F) + 2A{F{-2B + P +t)~rBF')f 

We see all the modes obey a dispersion relation cj2 — ABk^ . The appearance of the factor AB is due to the fact that 
t and r are coordinate time and distance. In terms of the physical time and distance, the dispersion relation says all 
the modes propagate at the speed of light, which is exactly the same as in the case of f{R) theories. Although there 
are no new contributions to the propagation speeds due to the Chern-Simons term, inclusion of the Chern-Simons 
yields a new correction to the propagation properties even when W' — 0. We find that the off-diagonal term ki^ is 
nonvanishing, 

_ 87r£(^ + ifMlYFs [rBA' [rF' + 2F) - 2A {rBF' + 2{B - l)F)) 
^^~ {21 + 1) (^3 -f 2^2 _ £ _ 2) A^F^Fx ■ ^ ' 

We also find that di3, 631 and 77113 are nonvanishing in general. This means that SF and q are coupled. On the other 
hand, there is no direct coupling between /3 and q. The coupling between SF and q means that we cannot consider 
the propagation of the odd and even modes separately as we can do in the case of f{R) theories, which is a clear 
difference from f{R) theories. This shows the potential usefulness of using the distinct nature of mode propagation 
in f{R, C) theories for putting constraints on f{R, C) models using observations of gravitational waves from compact 
astrophysical objects, when these become available in the future. 

We can also evaluate the mass for each eigenmode. However, since each matrix element of m^ is too lengthy to 
obtain analytic expressions for the mass eigenvalues, we will make an assumption that the background is very close 
to GR, i.e., F — 1, A ^ B = I — ^ and also expand the eigenvalues in £ = ^ (weak field approximation). Under 
these assumptions, three eigenvalues are given by 



1 £2(^_^_^)2^^_2(£2+^-2)^(^2_^^^2)r 6U{i^ + 3P + 7P + 9P-4:)rlF^ 



3Fa 3{P+e~ 2f r3 {P+i- 2f r^F^ 

i{e + i) e^e + ifrg 



O (e^) , (79) 



+ 0{e'), (80) 

+ 0(e^). (81) 



P {P + t- 2y r 

e^+i + i I92i {P + -it + 7P + 9P - 4) PgF^ 3 



p {p + e- 2y r^Fl 

Since ml is inversely proportional to Fx at leading order, this mode corresponds to the scalar graviton that exists in 
the general f{R) theories. To avoid the tachyonic mode, we need to impose a condition Fx > 0. The Chern-Simons 



^ There is also another reason the condition F > is required. As we can see from Eq. ^, GnF ^ gives the effective gravitational 
constant. Therefore, negative F yields repulsive gravitational force, which clearly contradicts our experiences. 
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corrections, Fg, appear in ml and rrig, but only in a combination with r^. This means those corrections are important 
only in the vicinity of the BH and are suppressed compared to the standard terms that exist in GR far from the BH. 
As we mentioned earlier, the results we derived do not apply to £ = and £ — 1 modes since we cannot solve ho 
and hi to go from Eq. (1701) to Eq. (|7ip . In the following, we consider the £ = and £ — 1 cases, separately. 

B. Monopole perturbation: 1 = 



The procedure we took from Eq. (|42|) to Eq. (j43l) cannot apply for the monopole case since bi identically vanishes 
for i = 0. Therefore, let us start again from Eq. (|^^ . For the monopole case, the odd modes and a identically vanish, 
which drastically simplifies the Lagrangian. In addition to this, the following relations among the background- 
dependent coefficients hold for £ — 0, 

fajb^y ajbs 0764 , . 



^2 / ^2 

Using these relations, we find that Eq. (|^^ reduces to 



C ^ Ho (a2H2 + ^SF + a^SF'] + —Hi ( a2i?2 + ^SF + ajSF'] + ciH^ + C3H26F' + C4H2SF 
\ 02 J aT \ 62 J 

+csH26'F + fi6F'^ + f26FSW + f36W^. (83) 

We can again eliminate 6W by using its equation of motion and the remaining fields are {Ho, Hi, H2,SF). For Ho 
and Hi, they appear only linearly. Therefore, they are Lagrange multipliers and their equations of motion yield 
constraints among the other fields. At first glance, we derived two constraints seem to over-constrain H2 and 6F. 
However, the two constraints are not independent of each other and we get only a single constraint which is given by 

a2H2 + ^-^SF + a7SF' ^~Ci, (84) 

02 

where Ci is an integration constant. This gives H2 in terms of 6F. 

It is important to notice here that the gauge conditions /3 — 0, K = Q and G = we have used to fix the even 
modes do not work for the £ = case since hta, h^a and the second term in Eq. (|26|) identically vanish. A condition 
that remains meaningful \s K — Q,hy which Rim is completely fixed and we still have another gauge transformation, 
Tfm = To- We can always set Hi = by properly choosing Tq- However, this condition still allows a remaining gauge 
degree of freedom of the form Tq = C2{t)A{r), where C2{t) is an arbitrary function. Correspondingly, Ho contains a 
gauge mode which is given by 2(^2 (i). Variation of Eq. ([83|) with respect to H2 gives a following equation: 

a2i?o = 2C1H2 + c^SF' + Ci5F. (85) 

Since H2 is written in terms of 5F by Eq. ([M|) . Hq is also determined by Eq. ((55|) once 5F is known. A homogeneous 
solution of a2i/o — is Ho — C2,{t), where C^it) is an arbitrary function, is a gauge mode and hence can be set to 
zero. 

Putting this constraint (|84|) into Eq. ([SS)) yields a Lagrangian which consists only of 5F, 



C 



vi U'f'^ - ABSF''^] + V2SF'^ + C1V3SF, (86) 



where vi , f 2 and U3 are given by 



2Att Mlr^F 

T2 ' 



VI = ^^. S _. , (87) 



/AB [rF' + 2Fy 
2nMj,r^VA vrMj, {rAA' {9rB' + iB) - 9r^BA'^ + 2A^ {3r^B" + 6rB' + 2B - 2)) 



V2 - 

VBI'x 2B'' 

2itMI Ua^ (srFB' {rF' + 2F) + 2B{F- rF'f) - rBA' {rF' + 2F) {rA'{rF' + 2F) + 2A{4rF' + 5F)))) 

^3 = ^ — ^ -^= 5 m 

A3/2y/B{rF' + 2Ff 

Since V2 has a very long expression, we have kept only terms which either diverge or remain finite in the GR limit 
Fa — )■ 0. The last term which is linear in SF acts as a source for SF. It is clear from the final Lagrangian that there 
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is only one propagating mode which does not exist in GR. This mode appears when F is allowed to fluctuate, as in 
the case of f{R) theories. The propagation speed is the velocity of light. The no-ghost condition vi > is satisfied 
if F > 0. Therefore, as is the case with i > 2, F > is the no-ghost condition. In order not to make the mode 
tachyonic, we need to impose another condition F\ > 0. 

Let us finally comment on the physical meaning of the other mode, which is accompanied by the constant Ci in 
Eq. ([84|) . This can be understood by considering the case of GR. In this case, we have SF — and 02 = —AMpirrA 
and ci = 2M|,7r. Then, Eq. 1^ becomes 






^^ = 4M|i^- (^«) 



Putting this into Eq. ((85)) gives 



III 



^" = 4Ml^ - ^- (^^) 



These solutions correspond to a metric. 



2Gn{M + C,)\ _,^2 , f, 2Gn{M + C,)\-\ 2 , 2, ^^2 , „• 2^ ^ 2^ 



Therefore, Ci just shifts a parameter of the background solution into another one. Physically, Ci represents shift of 
the BH mass. 

C. Dipole perturbation: £ — 1 

Let us start from Eq. ([68]). For the dipole case, we find that the relation q2qe — \<l\ = holds. This means only 
a linear combination of 5F and j3 can be a dynamical variable. To construct a dynamical variable, let us do a field 
transformation from 5F to "0 by 

i, = 5F+-^p. (93) 

2(72 

When written in terms of the new variables, no /3 appears in the Lagrangian. At the same time, we find all the /3' 
terms go away. This is due to another relation, q^^q-j — igf = 0. Now, since /3 is an auxiliary field, we can eliminate it 
by using its equation of motion. After this, the Lagrangian consists of tp, /ig and hi, which can be formally written as 

C = si (V'^ - ABii^'^^ + S2il>^ + - (s3 + — ) K + S3,h'^ + Sihiho - ^Sihih'^ + s^hj + sehl. (94) 

Interestingly, there is no coupling between the odd and even modes. 

Let us first focus on the even mode tp. We find that ip propagates at the speed of light. A coefficient in front of the 
kinetic term is given by 

2nMj,r^.f^{rBA' {rF' + 2F) - AA{B - l)Ff 

si = ^- o. (95) 

AFirBA' [rF' + 2F)-2A{rBF' + 2{B-l)F)f 

Therefore, as is the case with ^ > 2, F > is the no-ghost condition. The mass term is very complicated in general, 
but if the background is very close to that in GR, we find 

S2 = -^^. (96) 

In this case, the mode is not tachyonic if F\ > 0. 

Let us next consider the odd modes Hq and hi. The equations of motion for them are 

(2S3/10 - ^«4/ii) ' = S4A1 + - (4 + ^) ho, (97) 

f 

54/io — 7:^4^^o + '^sqHi — 255/11. (98) 
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It is important to notice that the Regge- Wheeler gauge /i2 = we have used to fix the gauge for the odd modes does 
not work for £ = 1 case. This is simply because hab for the odd modes identically vanishes ioi £ — 1 [see Eq. (^0)) ]. 
Therefore, there is still one gauge degree of freedom in the odd modes. We found it convenient to set hi — 0. However, 
as is clear from Eq. pop. the gauge is not completely fixed yet and we have a remaining gauge which is represented 
by A = Ci(t)r^, where Ci{t) is an arbitrary function. 

Under the gauge condition hi ~ 0, Eq. (jQT]) and Eq. (j98| become 



{ssh'.y ^1(4 + 7) ho. (99) 

^0 - ^K = 0. (100) 

The general solutions of Eq. (|100p are given by 

ho = C2{ty+yir), (101) 

where C2{t) and y{r) are arbitrary functions. Clearly, the first term C2(t)r^ is also a solution of Eq. ([M)) . But, this 
term is a gauge mode because it can be removed away by using the remaining gauge degree of freedom A — Ci{t)r^. 
On the other hand, the second term y(r) represents a physical solution. The form of y{r) is determined by the 
condition that y{r) must satisfy Eq. ([5^ which can be rewritten as 

y"+'-^y'--(- + '-^]y = 0. (102) 

S3 r \r S3J 

Since this is a second order differential equation, there are two independent solutions. Of course, one solution is a 
gauge mode which is proportional to r^ and can be absorbed into the first term of Eq. (|10ip . What we want is the 
other solution. To obtain a formal expression of the desired solution y, let us construct a following quantity out of y 
and r^: 

n = S3ir^y'^2ry). (103) 

By using Eq. (|102p . we can show that il' = 0. Setting il — const allows us to formally solve Eq. (|103|) in terms of y, 

where tq and Di are constants. Although there are two free constants we can choose, the shift of y due to shift of ro 
from one value to another only results in the change of the first term of Eq. ()101|) . which can be gauged away. On 
the other hand, solutions with different values of Di are not connected by the gauge transformation and hence Di 
represents a physical quantity. The explicit expression for S3 is given in the appendix [Cl 

Since y in Eq. ()104|) depends only on r and corresponds to the perturbation of the t — if metric component, 
this perturbation represents a stationary spacetime which slightly deviates from the static spacetime. Physically, it 
represents a spacetime metric around a slowly rotating BH, where the angular momentum of the rotating BH is taken 
into account as the perturbation. Therefore Di is related to the angular momentum of the BH. 

For demonstration, let us check that Eq. (|104p in the case of GR actually gives a Kerr metric expanded to first 
order in the angular momentum J, 

/ T \ cJt '2t T 

dsl^,, ~ - ( 1 - ^ ) rff2 + _ + ^2/^^2 ^ gjj^2 g^^2^ ff_ ^^^^2 q ^^^^^ /-^Q5^ 

V r / 1 Mr 

r 

where M is a mass of the BH. For the case of GR, we find .S3 = 4^Afp. Then, Eq. (|104p becomes y = — . ^U which. 



in terms of ht^, can be written as 



V = ^sin^^. (106) 



This actually coincides with the last term of Eq. (jlOSp for Di — —2 J. 
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V. CONCLUSION 

We have studied linear perturbations around the static, spherically-symmetric spacetime for general f{R,C) the- 
ories, where C is the parity violating Chern-Simons term. By explicitly constructing the second order action, we 
showed that one odd mode appears in the action as a quadratic in its second time derivative. Irrespective of its 
sign, this results in an Hamiltonian that is not bounded from below. Therefore, the static and spherically symmetric 
spacetime is unstable in general f{R,C) theories. This gives a strong limit on any phenomenological gravitational 
model which violates parity. 

We also showed that either R — const or qhqq = for the background metric is a necessary and sufficient condition 
to avoid the instability mentioned above. For such theories, the number of propagating modes for i? > 2 is three, one 
from the odd and the other two from the even. Unlike in the case of f{R) theories, those modes are coupled, which 
can be used as a distinctive feature to test the parity violating theories from observations. All the modes propagate at 
the speed of light. The no-ghost condition is ^ > and the no-tachyon condition is ^^ > 0, which are the same as 
in the case of f{R) theories. For the monopole perturbations {£ = 0), we showed that there is one propagating mode 
and one nonpropagating mode. Since the propagating mode amounts to the fluctuation of ^4, it does not exist in 
GR, but exists in f{R) gravity. It propagates with the velocity of light. On the other hand, the nonpropagating mode 
amounts to the shift of the mass of BH. This mode also exists in GR, as it should be. For the dipole perturbations 
{£ ~ 1), we found that the odd and the even modes completely decouple. The odd mode does not propagate and 
depends only on r. Physically, it corresponds to a slowly-rotating BH solution whose metric is linearized in its angular 
momentum. On the other hand, there is one even mode which does not exist in GR. This mode also propagates at 
the speed of light. The no-ghost condition and no-tachyon condition for £ = and £ ~ 1 are the same as for £ > 2. 
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2ttMI {rBA' {rF' + 2F) + A {ArBF' + F {2B + f + l))) 2^TM|,r^/AB [rF' + 2F) 

<^l = TT^ , r— — , 0-2 — 



as 



05 



{2(. + l)^/AB 
27r^(^ + l)M|,./4 [rFB' + 2B {rF' + F)) 



A7re{£+ 1)M^fVAB 
2er + r ' "" 2£+l ' 

27rM|, (rS^' (F - rF') + A {F {rB' + 2B + 2 {£^ + £ - l)) + 2rBF')) 



2^+1 



04 



(2£+1)FVAB 



27rM|,rJ4(7-B'+45) 
•^6 = .„ , , , 0.7 



A■nMlr'^^/AB 



as = 



21+1 ' 2^+1 

167r£(^ + 1)M|, (yl (VK' {2rB' + l^ + £ - 2) + 2rBW") - rBA'W) 

{2£+l)r^A ' 



8tt£{£ + l)Ml {rBA'W - 2A {rB'W + B {rW" - W'))) __ 167r£(£ + 1)M^BW' 



{2£+l)rA 



2^+1 



8Tr£{£ +l)Mj>{2A {rB'W + B{rW" + W'))- rBA'W) a^ 2a-j 

""- {2£+\)rA ' '^i^^-aio, 61-^, ^2-- — , 

, A' 2a2 , „, , 2aio , r ^^ ^ £ - 2 201261 

03 = -T^a?, 04 = — — , 65 = -26i, 66 = —-, br = --be, bs = ^ aio, 09 = , 



ci = 



C3 



ttMIJ^ {rA' {rF' + 2F) + 2A (2rF' + F)) 
2^+1 ^ 

27rM|,r./f (rA' + 4A) 



C2 = 



27r^(^ + 1)M|, Jf (rFA' + 2A (rF' + F)) 
2£r + r '' 



2^+1 



C4 



27rM|, (A (F (3rB' + 2B + 2 (^^ ^ ^ - l)) + 2rBF') - rBA' {rF' + F)) 



(2^+1)FVAB 



87r£(^ + 1)M|,S {rA' - 2A) W 
'''' ~ {2£+l)rA ' '^'^ 

2a4 flio 2aio 

rt A rA 



2c5 62 , 

C7 = -C5, Cg == — -, rfi 



r 
2aio 



2B' 



65 04 2a4 



F ' 



167r^(^ + l)MlBW' {2£{£ + 1)A - (^2 _^ ^ - 2) rA') 



, _ aio , _ ^(^+1) 
"9 — — T' ^10 — 5 aio, ei 



r2 ' " (2£+l)r3^ 

2ti£{£ + 1)M|, (A (F {rB' + 2B + e + £-2)+ 2rBF') - rBFA') 
(2<'+l)r2A3/2/B 



_ 46i _ 167r^(£ + 1)M2 (2A (rBF' + 2{B - 1)F) - rSA' {rF' + 2F)) 

^2-—. 63- ______ , 



64 



327r^(^ + 1)M|, (rSA' (r2F'2 + rFF' + F2) + ^ {-2r^BF''^ + r{l - 2B)FF' - 2{B - l)F^)) 

(2£+l)r3AF2 



„, , , {£-l)£{£+l){£ + 2)F 

65 = -26i, 66 = Oi, 67 = 61, eg == —^ 07, eg = 63, 



27rM|,r2VK^y'4 
(2^+1)(F2-FaM/«)' 



2r4 
2F5 f — ^'^ f 



(Al) 
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Appendix B: Expressions of pi, p2, • • • • 

Since most pi,- ■ ■ have very long expressions, we assume that F', W', F\, Fs and Ws are small quantities and 
truncate at the leading order, pi, p2 a-nd p^ are exact. 



2 TI//2 



Pi 



Pi 



32Tre{i + l)Mf,W 



(2£+l)F(4) 



A\3/2 



P2 



2pi 



P3 =Pl, 



32ne{i + \)MlW' 



P5 = 

P6 = 

P7 ^ 

P9 = 

PlO = 

Pl2 = 

Pl3 = 

Pl5 = 

Pie = 

Pl9 = 

P20 = 



(rA^BA' (W {-3i{e + 1) (-2rB' + ^^ _^ £) 

(2^ + Ijr^A^F^^ {rBA' + A {-2B + P+ £)f 

-AB {3rB' - 8) - 205^) + 8rB {-2B + f + l) W") + r^AB^A'^ {W' {3rB' + UB - U{1 + 1)) + ArBW") 
-ir^B^A'^W' + A^ {;di^{(. + ifrB'W' + 45^ {W' {3rB' + 2{f +£-4)) - Mil + l)rW") 
-2£{e + l)B {W' {6rB' + 71^ + 71- 16) - 2l{£ + l)rW") + 8B^ {2rW" + W'))) , 
2TTi{e + l)MlF 

(2£+l)J| 



AirMJ^r^-'^ 



2d2 4/2 



{{2£^ + 21- l)r^B^A 



£ (2£2 + 3^+1) AF {rBA' + A (-25 + £2 + ^))" 

+2£{£ + l)rABA' (-3B + 2£^ + 2l-l)+ £{£ + 1)^^ (_6^(^ + 1)B + 65^+^ (2^^ + A£^ + £ - l))) , 



87r^(^ + 1) (^2 _^ ^ _ 2) MlA^F 



Ps 



327r (^2 + ^ _ 2) M^W {rBA' + £{£ + l)Af 



{2£+l){Af^{rBA' + A{-2B + £^ + e)f' '" {2£ +l)AF {rBA' + A{-2B + £^ + £)y 

MttI{£ + 1) (^2 _^ ^ _ 2) MlBW {rBA' + £{£ + 1)A) 

{2£r + r) {rBA' + A {~2B + £^ + £)f 
Stt {£^+£- 2) MlrBJ^ {rBA' + £{£ + I) A) 



{2£ + 1) {rBA' + A {-2B + £^ + £)Y 



Pii 



27r^(^ + l)Ml,F 



AirM^r^BJ^ 



£ (2£2 + 3^ + 1) i^ {rBA' + A (-25 + £^ + £))' 



(2^+1) 
((2£2_^2^-l)r2B2^'2 



+2£{£ + l)rABA' (-3B + 2£^ + 2^ - l) + £{£ + 1)^^ {-6£{£ + 1)5 + 65^+^ {2£^ + 4^^ + ^ - l))) , 



Pu = 



8ti£{£ + l){£'^+£- 2) NPpA^F 
{21 + 1) (4)'/' {rBA' + A (-25 + £^ + £)f 
Mt:£{£ +\){£'^+£- 2) Mj,AB'^W' {rBA' + £{£ + I) A) 
{2£r + r) (r5yl' + yl (-25 + £2 + £))2 



327r (£2 + ^ - 2) MleBW {rBA' + £(^ + 1)^)^ 
{2£ + 1)F (r5A' + A (-25 + ^2 _^ £))2 



Stt (£2 + ^ - 2) M|,rA3 (r5A' + £{£ + l)yl) 
(2£ + 1) {^f^ {rBA' + A (-25 + £^ + £))" ' 

327r£(£ + l)Mf>Fs {rBA' - 2A(5 - 1)) 
(2£+l)r2AFA ' 



Pl7 



47r£(£ + 1)M|,F 
(2i?+l)7| 



Pis 



8it£{£+1)MIF 
{2£r + r)Jl ' 



327r£(£+l)M2 5 



{rABA'W' {{£^ + £ - 2) rB' + A£{£ +l)B + £ {£^ + 2£^ - 5£ - 6)) 



{2£ + l)r2 {rBA' + A (-25 + £^ + £))' 

+ {£^+£-A) r'^B'^A'^W' + A^ {W' {£{£ + 1) ((^^ + ^ _ 2) r5' - 2£{£ + 1)) 

+25 (-2 (£2 + ^ _ 2) rB' + 1'^ + 2£^ + b£^ + A£ - A) - A£{£ + 1)52) - 4 (^2 ^ ^ _ 2) rB^W")) 

32n£{£ + l)Mj,Fs {2A{B - 1) - rBA') _ Mn£{£ + 1)M|,5W^' 

P21 — ^0/7 I i\ 9 ■< 7^ ' P22 — 



{2£+l)r^AFx 



P23 = 



{2£ + 1) (r5^' + A (-25 + £2 + ^)) 



{2£ + l)r2 
2 (rA5^' ((^2 + ^ + 1) r5' + (4^2 + 4£ - 2) 5 - £ (2^3 + 4^2 _^ ^ _ 1)) 



+ (-2^2 -2^+1) r^B^A'^ +A'^(b({£^+£-2)-2{£^+£+1) rB' 
+e {e^ + 2^2 + 2£ + 1) rB' + 2 (£2 + £ _ 2) 52)) , p24 



27r£2(£+ 1)2^2 



(2^+1) 



(Bl) 
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6ATr£{e + l)Mj,F,F' 
^^^ ~ {2l+l)r^Fx ' 

647r^(£ + 1)M|, ((^4 + 2£^ +£^-A)W' ~ 2r (rVK^^) „ 2W")) 27r£ (^^ + 2^^ _ ^ _ 2) Mj, 

P26 = ^o«.^ , o/,9 — ^7 — 7;n^ : P27 - 



(2^3 + 3£2 _ 3^ _ 2) r3 ' """ (2£+l)r2 



27rM?,r2 
P2S ^ 



(2^ + 1)Fa ' 
P29 = ^ — :.(6(f+£+l)A^-4(£^ + e+l)A^ 

+ (^^ + 3^^ + 6^'' + 71^ + 7£^ + M+2)A- 3^1* - (^^ + £ + l)^) , 
87r£ (£3 + 2^2 _ ^ _ 2) M|,A (-3 (^^ + £ + 1) ^2 _^ ^^2 ^ ^ ^ ^^^2 ^ + 3^13 _ ^2 _ ^ _ ;l) 
^^° ~ {2£+iy{-SA + £^+£+lf ■ 

Appendix C: Expressions of S3. 

S3 is given by 

zi 
S3 = — , 

Z2 

where Zi and Z2 are given by 
zi = 47rAf|, ('64r3A2F2F' (rSA' (rF' + 2F) - 2^ (rSi^' + 2{B - l)F)f 

— ^ {-r^AFxA' {3F^B' + 2BFF' + ArBF'^) + 3r^BF^FxA'^ - 2A^ {F^ (3Fa {r^B" -2B + 2) + r^F') 

„2 D I?. I?'2 o„/' D 1 \ I? I?. I?'^^ / '}0^2 o2 ^/2 /„ t;i/ , o I?^2 /' ei2 



-ir^BFxF'^ - 2r{B - l)FFxF')) {-■i2r^B'^A'^ {rF' + 2F) {F^ - FxWs) 

+128rABA' {rF' + 2F) {rBF' + 2{B - l)F) {F^ - FxWs) + A^ {-128r^B^F'^ {F^ - FxWs) 
-512r{B - l)BFF' {F^ - FxWs) - 512{B - ifF'^ (F^ - FxWs) + r'^F^Fx))) , 

[a 

Z2 = Sr^A^FxJ — {r^AFxA' {3F^B' + 2BFF' + ArBF'^) - ir^BF^FxA'"^ 

+2A^ {F^ {3Fx {r^B" - 2B + 2) + r^F') - 4r^BFxF'^ - 2r{B - l)FFxF')) . 
In particular, for f{R) theories in whicli Fs — Ws — 0, S3 is drastically simplified to 

.3 - '-^. (CI) 
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